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Abstract 

We apply a technique, due to Stephani, for generating solutions of the 
Einstein-perfect-fluid equations. This technique is similar to the vacuum so- 
lution generating techniques of Ehlers, Harrison, Geroch and others. We start 
with a "seed" solution of the Einstein-perfect-fluid equations with a Killing 
vector. The seed solution must either have (i) a spacelike Killing vector and 
equation of state P = p or (ii) a timelike Killing vector and equation of state 
p + 3P = 0. The new solution generated by this technique then has the same 
Killing vector and the same equation of state. We choose several simple seed 
solutions with these equations of state and where the Killing vector has no 
twist. The new solutions are twisting versions of the seed solutions. 



I. INTRODUCTION 



In solving Einstein's equation it is usual to treat spacetimes that have one or more Killing 
vectors. The existence of a Killing vector reduces Einstein's equation from a 4 dimensional 
set of equations to a 3 dimensional set. Instead of the spacetime metric one considers an 
equivalent set of variables: the norm A and twist potential u of the Killing vector and a metric 
jab on the 3 dimensional manifold of trajectories of the Killing field. The vacuum Einstein 
equations then become a set of equations for the quantities (A, u, ja,)- Remarkably, as shown 
by Ehlers [JTJ , Harrison M , Geroch pj and others HI] these equations have a symmetry that 
allows us to generate new solutions: given (A, u, 7 a &) solving the equations then (A, {jj, 7^ 
is also a solution where (% uij is a simple algebraic function of (A, u). These results have 
been generalized to the Einstein-Maxwell equations @,|5l,|6j and have been applied to find 
many new solutions. (The general invariance transformation was first published in |4| (in 
the vacuum case with two Killing vectors) and in || (in the Einstein- Maxwell case with one 
Killing vector). Ehlers [|IJ found a special case of the general vacuum transformation while 
Harrison found a special case of the general Einstein-Maxwell transformation.) 

The vacuum solution generating technique was generalized to the case of perfect fluids 
by Stephani [[7]. This technique is very similar to that for the vacuum Einstein equation. 
However, it turns out that not all equations of state are suitable for this sort of solution 
generating. Only two equations of state are compatible with this technique: P = p for a 
spacelike Killing field and p + 3P = for a timelike Killing field. 

In section 2 we define our notation. Sections 3 and 4 contain the new solutions that we 
have generated using this technique with section 3 treating the P = p case and section 4 
treating the p + 3P = case. Section 5 contains a brief discussion of our results. 
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II. FORMALISM 



Our notation is essentially that of refs. [[|[7|]- We use signature (—,+,+,+) and units 
where c — G — 1. Let (M,g a iy) be a solution of the Einstein-perfect fluid equations with 
energy density p and pressure P. Assume that g a b has a Killing vector £ a . Define the norm 
A and twist uo a of £ a by A := £ a £ a and 

oo a :=e abcd eVX d • (2.1) 

We assume that either (i) the Killing vector is orthogonal to the fluid four-velocity and 
P = p or (ii) the Killing vector is parallel to the fluid four-velocity and p + 3P = 0. From 
this assumption it follows that R a b£, b = 0. It then follows that there is a scalar uj such that 
w a = V tt w and that there are forms a a and (3 a satisfying 

V [o a 6] = -e o6cd V c £ d , (2.2) 

i a a a = u , (2.3) 

V [a/ 6] =2AV a & + uje ahcd WX d , (2-4) 

CPa = to 2 + A 2 - 1 . (2.5) 
For any constant 9 define A and r\ a by 

A/A := cos 2 fl + (uj 2 + A 2 ) sin 2 fl - 2uj sin 9 cos 9 , (2.6) 

•q a :=\~ 1 ^ a + 2 cos 9 sin 9 a a - sin 2 ^^ . (2.7) 
Then the new metric is given by 

9ab=^ (gab ~ ^LZb) + >^VaVb ■ (2.8) 

A 

This new metric is also a solution of the Einstein-perfect fluid equations with the same 
equation of state and the same Killing vector. 

One of the uses of solution generating techniques is that they can produce twisting 
solutions from seed solutions that are twist free. We now consider how the process for 
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finding the new metric simplifies in the case where the seed metric is twist free. In this case 

u = and V a £ b = A" 1 £ [6 V a] A. Define F := cos 2 6 + A 2 sin 2 #. Then A = X/F. The equation 
for p a can be solved as p a = (A — A -1 ) £ a . The equation for a a simplifies to 

V [a a b] = e abcd i d V c \ (2.9) 

and ^ a a a = 0. The one-form r] a is then given by 

Va = A" 1 ^ + 2 cos6 sm6a a . (2.10) 

The new metric is then given by 

9ab = F (g ab - X^iaib) + ^VaVb ■ (2.11) 

From now on we will assume that the seed solution is twist-free. 



III. P = p FLUIDS 
A. Fluid variables 

In describing the new solutions we will concentrate on two properties of interest: singu- 
larities and fluid flow. The solution generating procedure often results in singularities on 
the axis. For each of our solutions we will consider whether such a singularity exists. A 
perfect fluid is described by the energy density p, pressure P and fluid four velocity u a . The 
covariant derivative of the fluid four velocity is decomposed as 

VfcWa = -a a u b + to ab + a ab + 6 (g ab + u a u b ) /3 . (3.1) 

Here a a := u b V b u a is called the acceleration. The tensors a ab (called the shear) and u ab 
(called the vorticity) are orthogonal to u a with u ab antisymmetric and a ab symmetric and 
trace-free. The scalar G is called the expansion. Under solution generating with the P = p 
equation of state the fluid parameters transform as follows: 

p = p/F , (3.2) 
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u a = VFu a , (3.3) 
® = -jf ( e + i ? ~ 1 Asin 2 0u o V o A) , (3.4) 



a a = a a + F" 1 Xsin 2 9 (V a A + u a u b V b XJ , (3.5) 

uj ab = \/F io ah . (3.6) 

There does not seem to be a simple transformation rule for a a b. 

A fluid with the equation of state P = p is called "stiff matter." In such a fluid the 
adiabatic sound speed is equal to the speed of light. If the fluid is irrotational then the field 
equations are the same as those of a minimally coupled scalar field with timelike gradient. 



There is a technique for generating stiff matter solutions from vacuum solutions. [10 



B. Robertson- Walker seed metric 

Consider the seed metric 

ds * = _ df 2 + t 2/3 ^2 + dr 2 + r 2^2] (3 7) 

This is a flat Friedmann-Robertson- Walker metric with P = p. The fluid four velocity is 
u a = —V a t and the energy density is p = 1/ (247rt 2 ). The shear, vorticity and acceleration 
all vanish. The expansion is = 1/t. The metric is type O. This metric has more than one 
Killing vector. We first consider the case = (d/dz) a . Then we have A = t 2/3 and uo = 0. 
It then follows that F = cos 2 8 + t 4//3 sin 2 ^. The equation for a a is then 
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V [a a fc] = - - rV [a rV 6 ]0 (3.8) 
(as well as ^ a a a = 0). A solution is 

a a = -\r 2 V a (i> . (3.9) 

The one-form r] a is then given by 

2 

Va = V tt z — - cos 9 sin 9 r V a • (3.10) 

o 



The new metric g a b is then given by 

d ~ s l =F [_ d f2 + f2/3 ^ r 2 + r 2 #2 y 

This new metric has no singularities on the axis. The acceleration and vorticity vanish. The 



+2/3 , 2 \ 2 

+ — f cfe - - cos# sin#r 2 d0J . (3.11) 



expansion is 



e= * (V 1 + -F-H^hi^e] . 

y/F V 3 / 

The shear has components in the rr, 00, zz and £0 directions. Its square is 



/32X t 2 / 3 sin 4 fl 
cr a ab = — I 



(3.12) 



(3.13) 



V277 F 3 
The metric is type I. 

We now use the same seed metric but with the Killing vector £ a = (<9/<90) a . Then we 
have A = t 2/3 r 2 and u = 0. It then follows that F = cos 2 9 + t 4/3 r 4 sin 2 #. The equation 
for a a is then 



V [a a b] = -rV [a rV h] z + 2t 1/3 V [a t V b] z 



(as well as £ a a a = 0). A solution is 



1 2 3 



,3 2 
The one-form r) a is then given by 

r?a = V a + (Jr 2 + 3t 4/3 ^) cos9 sin9V a z 

The new metric g ab is then given by 



(3.14) 



(3.15) 



(3.16) 



+ 



ds 



£2/3 r .2 



7,2 



+ 



- dt 2 + t 2 ' 3 (dz 2 + dr 2 ) 
2 



~ 2 + 3 t 4 / 3 



cos# sin 9 dz 



(3.17) 



This metric has no singularities on the axis. (There is a conical singularity if the range of 
is 2n. However, this singularity can be removed by choosing the appropriate range of 0). 
The vorticity vanishes. The expansion and acceleration are given by 



9 = ^= (V 1 + - F-y/ 3 r 4 sin 2 
y/F V 3 



(3.18) 



d a = 2F- 1 t 4/3 r 3 sin 2 # V a r . (3.19) 

The shear has components in the zz, rr, <jxj) and zip directions. Its square is 

8 \ r 2 t 2 / 3 sin 2 



a a n 



27) F 3 



27cos^ + 4r t W0 



(3.20) 



The metric is type I. 



C. Tabensky-Taub seed metric 

We now consider a different seed metric: The Tabensky-Taub metric |§ 

ds 2 = -Vdt 2 + Vdz 2 + z(dx 2 + dy 2 ) (3.21) 

where 

V = z- 1/2 exp(^a 2 z 2 ) . (3.22) 

This is a solution of the Einstein-perfect-fluid equations with 

a 2 



The fluid four-velocity is u a = — VV V a t. The expansion, vorticity and shear vanish. The 
acceleration is 

1 dV . 

a ' = 2v7I v " z ■ (3 ' 24) 

The metric is type D. 

First we consider the Killing vector £ a = (d/dx) a . Then we have A = z and thus 
F = cos 2 9 + sin 2 6. The equation for a a is V[ a afe] = — V[ a y V&]£ (as well as (, a a a = 0). 
The solution is a a = — y V a t. It then follows that 

f]a = V a £ — 2 cos9 sin 9 yV a t ■ (3.25) 

The new metric is then given by 

ds 2 = F (-Vdt 2 + Vdz 2 + zdy 2 ) + J (dx - 2 cos 9 sin9ydt) 2 . (3.26) 
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Both the seed metric and the new metric have singularities at z — 0. However, the new so- 
lution does not appear to have any additional singularities. The expansion and the vorticity 
vanish. The acceleration is 

/ 1 dV zsm 2 6\ „ 

[zv 17 + —) ■ (3 ' 27) 

The shear has components in the xy and ty directions. Its square is 

nh 2cos 2 #sin 2 # 

° ab ^= FW ■ (3-28) 

The metric is type I. 

We now consider another Killing vector of the Tabensky-Taub metric. Using polar co- 
ordinates for the x — y plane we find that the seed metric is 

ds 2 = -Vdt 2 + Vdz 2 + z(dr 2 + r 2 d<f) 2 ) . (3.29) 

Using the Killing vector £ a = (<9/<90) a we find A = zr 2 and u — 0. Therefore F = 
cos 2 9 + z 2 r 4 sin 2 ^ . The equation for a a is 

V [a a b] =rV [a rV b] t + 2VV [a tV b] z (3.30) 

(as well as £, a a a = 0). A solution is 

a a = 2VtV a z - rtV a r . (3.31) 

The vector i] a is then given by 

r] a = V a (f) + 2 cos 6 sin 6 {2VtV a z - rtV a r) . (3.32) 
The new metric is then given by 

ds 2 = F [-Vdt 2 + Vdz 2 + zdr 2 ) 

2 

+ ^- {d(j> + 2cos# sin^ [2Vtdz - rtdr]) 2 . (3.33) 
F 
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This metric does not have any singularity on the axis (though one must choose the range of 
<p to avoid a conical singularity). The expansion and vorticity vanish. The acceleration is 



/ 1 dV zr i sm 2 9\ „ 2z 2 r 3 sm 2 6 „ 

The shear has components in the rr, r(f), rz, zz and <fiz directions. Its square is 

2 r 2 cos 2 # sin 2 fl (r 2 + AzV) 



F 3 V 



The metric is type I. 



D. Allnutt seed metric 



(3.34) 



(3.35) 



Our next seed metric is due to Allnutt. It is 

dt 2 



ds 2 = z 2 



,,2 t 2 (l+t 2 ) 



dy* 



+ dz 2 



(3.36) 



where the function / is / = t 2 ^ (1 + t 2 ) 1 ^ and (5 is a constant. This spacetime is a solution 
of the Einstein-perfect fluid equations with 

/3(l-/3) 



P = p 



(3.37) 



8nz 2 t 2 (1 + t 2 ) 

and u a = — z (1 + 1 2 ) 1 ^ 2 V a t. The vorticity vanishes. The acceleration is a a = z~ 1 'V a z and 
the expansion is 



e 



2t 2 + 1 



ztVTTP 

The shear has components in the xx, yy and zz directions. Its square is 



(3.38) 



3z 2 t 2 (1 + 1 2 ) 



t 4 + t 2 + 1 + 3/3(/3 - 1) 



(3.39) 



The metric is type D. 

Using the Killing vector £ a = (d/dx) a we find A = z 2 f and thus F = cos 2 9 + z 4 f 2 sin 2 #. 
The equation for a a is 



W [a a b] =2z 2 tV [a tV b] y + 2z (t 2 + (3) V [a zV b] y 
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(3.40) 



(as well as £ a a a = 0). A solution is ot a = z 2 (t 2 + (3) V a y. It then follows that 



Va = V a x + 2 cos sin z 2 (t 2 + (3\ V a y . 
The new metric is then given by 
d~s 2 = F 

+ -TT- 



(3.41) 



z 2 t 2 (l + t 2 



dt 2 + 



dx + 2 cos sin z 2 (t 2 + p) dy 



dy 2 + dz 1 



(3.42) 



Both the seed metric and the new metric are singular at z = and at t = 0. However, 
the new metric does not appear to have any additional singularities. The vorticity vanishes. 
The expansion is 

~ fl + ¥ ( 2t 2 + l z 3 ff sm 2 6 \ 

= i— [zwt¥) + -^—) (3 - 43) 

where an overdot denotes derivative with respect to t. The acceleration is 



i a = [z- 1 + 2F~ 1 z z fsin 2 e) V 



(3.44) 



The shear has components in the xx, yy, zz and xy directions. Its square is 

2 



<? CTab 



3z 2 t 2 (l+t 2 )F 



t 4 + t 2 + 1 + 3/3(/3 - 1) 



8 • 2n 2 P 

- sm 0z — - 
3 F 2 



t 2 + (3 (F_ 1 - 3/3 - t 2 \ 3cos 2 g 
~T~ \F + t(l + t 2 ) J + F 



(3-45) 



The metric is type I. 



E. Kramer seed metric 



The final P = p seed metric that we shall consider is due to Kramer. [0 It is 



ds 2 = e a2r2 (-dt 2 + dr 2 ) + r 2 d^ 2 + dz 2 



(3.46) 



where a is a constant. This is a solution of the Einstein-perfect fluid equations with 



P = p = — e 

87T 



2 



(3.47) 
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The four velocity is u a = — e° 2 ' r2//2 V t. The expansion, shear and vorticity vanish. The 
acceleration is a a = a 2 rV a r. The metric is type D. We use the Killing vector £ a = (d/d<p) a . 
Then we have A = r 2 and thus F = cos 2 9 + r 4 sin 2 6 ) . The equation for a a is Vua&i = 
— 2V[ a zVh]t (as well as ^ a a a = 0). A solution is a a = — 2zV a t. The new metric is then 
given by 



d~s 2 = F ' ' ' ' 1 ' ' 



2 



e a2r2 (-dt 2 + dr 2 ) + dz 2 } + — {d(f> - 4 cos 9 sin 9 zdtf . (3.48) 

The new metric has no singularities on the axis, though one must choose the range of <p to 
avoid a conical singularity. The expansion and vorticity vanish. The acceleration is 

a a = (a 2 r + 2F" 1 r 3 sin 2 ^) V a r . (3.49) 

The shear has components in the zip and zt directions. Its square is 



2 „-o 2 r 2 



i r e 



cos^sin 2 ^ 



a ab ~a ab = ^^ ^ — - • (3.50) 



The metric is type I. 



IV. p + 3P = FLUIDS 



A. Fluid variables 



We now consider the solutions with the equation of state p + 3P = 0. This equation 



of state occurs in the Einstein static universe in one of the limiting solutions of the 



Wahlquist metric |TJ| due to Vaidya RfOJ and has been suggested as the equation of state for 



a chaotic array of strings |TJ| . In generating solutions with the equation of state p + 3P 



we must choose the fluid four-velocity of the seed metric to be parallel to the Killing 
vector. The fluid four velocity of the new metric will then also be parallel to the Killing 
vector. This results in considerable simplification of the fluid parameters. The expansion 
and shear of both the seed solution and the new solution vanish. The density, four-velocity 
and acceleration transform as follows: 

p = p/F , (4.1) 
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u a = VFu a , (4.2) 

a a = (2F~ l cos 2 9 - l) a a . (4.3) 

Since we choose seed solutions with zero twist, it follows that the fluid vorticity of the seed 
solution vanishes. The vorticity of the new metric is 

u ab = 2 cos 9 sin 9 y/- X/F V [a a b] . (4.4) 

It then follows that the vorticity vector Q a := e a bcd ; u b £j cd is parallel to the acceleration vector 
a a . 



B. Vaidya seed metric 



Our first seed metric is due to Vaidya ]T5| and is given by 

d s 2 = -fdt 2 + - — dr ox - + r 2 (dip 2 + sm 2 i/jd(f) 2 ) . (4.5) 
(1 + cr 1 ) j v ' 

Here 

f . = 1 _ 2m v / 1 - cr2 (4 6) 

r 

and c and m are constants. This two parameter family of solutions contains the 
Schwarzschild metric and the Einstein static universe as special cases. The energy den- 
sity is p = — 3c//(87r) and the fluid four velocity is u a = —y/fV a t- The acceleration 
is 

«. = £ § VT ■ (47) 

The metric is type D. The timelike Killing vector is £ a = (d/dt) a . It then follows that 
A = — / and therefore that F = cos 2 9 + f 2 sm 2 9. The equation for a a is V[ a a b -\ = 
— 2m sin ip V[ a ip Vb}4> (as well as (, a a a = 0). A solution is a a = 2m cos^ V a </>. We then 
find 

Va = V a £ + 4m cos 9 sin 9 cos ip V a . (4.8) 
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The new metric is 



ds 2 = F 



dr 2 



;i + cr 2 ) / 



+ r z 



+ sin 2, d(f)' 



f 2 
— — (dt + 4m cos 9 sin 9 cos ipdcf)) 



This metric has a singularity on the axis. The acceleration is 

5.= (2F-W0- l) Yffr Var ■ 
The nonvanishing components of the vorticity are 

[J 

tiip<t> — ~ &H> — ~ 2m cos 9 sin 9 \ — sin ip 



(4.9) 



(4.10) 



(4.11) 



After the transformation the pressure and energy density are scaled by F _1 from their static 
values. As in the static case, on the horizon the pressure and scalar curvature are zero. The 
generated vorticity tensor vanishes on the horizon. The metric is type I. 

C. Kramer seed metric 



The next seed metric that we will consider is essentially due to Kramer. It is 



ds 2 



e 2x df + 



1 



e 2x sinh bx 



■- <K.r + e v "~ dz + — 



dx 2 



(4.12) 



sinh bx 

Here b is a constant. The energy density is p = 3b 2 e 2x /(327r) and the fluid four velocity is 
u a = — e x V a t. The acceleration is a a = V a x. The metric is type I. The Killing vector is 
£ a = (d/dt) a and therefore we have A = — e 2x . It then follows that F = cos 2 9 + e 4x sin 2 #. 
The equation for a a is V[ a afe] = 2V[ a 2 Vb}4> (as well as i a a a = 0). A solution is a a = 2zV a 0. 
We then find 

Va = V a t + 4 cos 9 sin 9z V a <p . (4. 13) 



The new metric is then given by 
F 



ds 



7:2 



e 2x sinh for 



^ x d6 2 + e-^ Ax dz 2 + 



dx 2 



sinh bx 



•2x 



6 2 

— (dt + 4 cos 9 sin 9 z dtp) 
F 

(4.14) 
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The solution generating process does not seem to have added any new singularities. The 
acceleration is 

a a = (2F- l cos 2 6 - l) V a x . (4.15) 
The nonvanishing components of the vorticity are 

Co Z(j> = - Cb^ = 2 cos 9 sin 6e x F~ 1 . (4.16) 

The metric is type I. 

V. DISCUSSION 



Exact solutions to the Einstein field equations are valuable as a tool in the exploration 
of behavior ranges allowed by the field equations. With the marked increase in numerical 
calculations, exact solutions are also useful as comparisons for approximate or numerical so- 
lutions and as checks for the development of computer codes. In general Einstein's equations 
for static spaces are far easier to solve than those for stationary spacetimes. One advantage 
of solution generating methods is that they can be used to find stationary solutions using 
simple static seed solutions. In this paper, we have used the method of Stephani to 
generate stationary solutions of the Einstein-perfect fluid field equation from simple static 
seeds. A possible disadvantage of solution generating is that often extra singularities, espe- 
cially on an axis, can be introduced. Some of the solutions that we generated have no new 
singularities and in some cases they appear to be non-singular. 

While solution generating techniques have been applied extensively to vacuum spacetimes 
and to solutions of the Einstein-Maxwell equations |§; with some notable exceptions PJlO 



there is comparatively little work on solution generating involving perfect fluids. The method 
of reference [JTO]] generates P = p perfect fluid solutions from vacuum solutions. It can be 



applied only to solutions with two spacelike Killing vectors and generates only solutions 
with zero vorticity. In contrast the Stephani method requires only one Killing vector and is 
compatible with nonzero vorticity. 
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There has been recent numerical work on vacuum cosmologies with two spacelike Killing 



vectors. W7\ In this work the known vacuum exact solutions provided useful test cases. The 



numerical work of reference JT7| could easily be generalized to the case of inhomogeneous 
perfect fluid cosmologies. The exact solutions that we have generated in this paper could 
be used as test cases for such numerical work. 



All of the seed metrics that we have used belong to the Perjes Class II ||18|| . Since 
this classification is preserved under the transformation [|l9j the generated metrics are also 
members of this class. For the P = — p/3 solutions, the generated solutions are all shear and 
expansion free and all have a non-zero vorticity that lies along the static acceleration. This 
is a characteristic of stationary Perjes II solutions. No assumption about algebraic speciality 
is made. In the P = p case, the fluid parameters of the new solutions are not constrained 
to be shear free or rigidly rotating. The action of the transformation in this case is on 
the spacelike Killing vector and all of the generated solutions as well as their seeds have 
zero vorticity while their shear, expansion and acceleration structure can be quite complex. 
Solution generating allows access to exact solutions that are difficult to find using other 
methods. These new solutions illustrate various behaviors allowed by the Einstein-perfect 
fluid equations and may additionally serve as aids in further numerical work. In some cases 
they provide insights into the static origins of existing solutions obtained through more 
arduous methods. 
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